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^ : Abstract 

: We report 75As nuclear magnetic resonance studies on a very clean hole " doped single - crystal 

*^ ' Bao 68^0 32Fe2As2 (T c = 38.5 K). The spin-lattice relaxation rate 1/Ti shows an exponential de- 

G ■ 

crease below T ~ 0.45T C down to T ~ 0.11T C , which indicates a fully-opened energy gap. From the 

o ; 

ratio (Ti) c /(Ti) a , where a and c denote the crystal directions, we find that the antiferromagnetic 
spin fluctuation is anisotropic in the spin space above T c . The anisotropy decreases below T c and 
' disappears at T — > 0. We argue that the anisotropy stems from spin-orbit coupling whose effect 

^t- : 

vanishes when spin-singlet electron pairs form with a nodeless gap. 
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The discovery of superconducting transition in electron-doped iron-arsenide 
LaFeAsOi-^Fa; provides a new route to high temperature superconductivity Re- 
markably, many other .RFeAsOi-aT-j, (R: rare earth) were synthesized and T c was raised to 
55 K in SmFeAsOi-ajFj. j^j, which is the highest among materials except cuprates. Soon 
after these works, the hole-doped BaFe 2 As 2 was also found to be superconducting |3]. The 
large single crystals of Bai_ x K x Fe2As2 are easy to obtain, which makes them a good system 
for studying many physical quantities. 

One of the most outstanding issues for a new superconductor is the symmetry of the 
electron pairs which is directly related to the paring mechanism. Nuclear magnetic resonance 
(NMR) experiments found the electronpairs to be in the spin-singlet state [4| and indicated 
the existence of multiple energy gaps 0,0]. The multiple-gap property is likely associated 
with the multiple electronic bands. The Fermi surfaces consist of two hole-pockets centered 
at the T point and two electron pockets around the M point {(]]. However, whether there 
are nodes in the gap function or not is still under hot debate. Angle-resolved photoemission 
spectroscopy (ARPES) jz| suggested fully opened gaps, but thermal conductivity [s[ Q 



measurements suggested nodal gaps, 
groups have led to opposite conclusions 



he penetration depth measurements by different 
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Theoretical 
ing candidate 



y, the sign-reversing s ± -wave model has been considered as the most promis- 



12l4l4 L but d wave or s wave with zero gap, and even a conventional s ++ wave 



was also proposed 15rll7|. It has been shown that the s ± -wave or a multiple-gap d-wave 
model can fit quite well the spin-lattice relaxation rate, 1/Tj , which shows a rapid decrease 

-l20j]. However, an important feature 



ay, 
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below T c with a hump structure at T ~ T c /2 
that 1/Ti should decrease as an exponential function of T expected for the s-wave gaps 
has not been observed so far, because of impurity scattering in the samples. The impurity 
scattering; can also alter other physical properties Q. Thus, the conclusions on the gap 
symmetry drawn so far are still controversial. Measurements in sufficiently clean samples 
are highly needed to resolve the issue. 

Here we report 75 As NMR study on a very clean single crystal Bao.6sKo.32Fe2As2 with 
T c = 38.5 K that is the highest among reports for this family. We obtained two pieces 
of evidence for fully-opened gaps. First, we observe an exponential decay of 1/T\ below 
T ~ 0.45T C down to T ~ 0.11T C . For the second piece of evidence, we find that the 
antiferromagnetic (AF) spin fluctuation (SF) is anisotropic in the spin space above T c , but 
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FIG. 1: (Color online) (a) 75 As-NMR spectra at a frequency of u /2ir = 55.1 MHz and T=100 K. 
The vertical axis for H \\ a is offset for clarity, (b) The T dependence of the Knight shift with 
H || a axis and H \\ c axis, respectively. The arrow indicates T c for H \\ a. 



the anisotropy decreases below T c and disappears at T — )■ 0. We argue that the anisotropy 
is due to spin-orbit coupling, whose effect vanishes at T — > because the electron pairs are 
in the spin-singlet state with nodeless gap. 

The single crystal of Bao.6sKo.32Fe2As2 was grown by using the self-flux method and 
characterized as discussed elsewhere [221 ] . Both dc susceptibility measured by a super- 
conducting quantum interference device and ac susceptibility measured by the NMR coil 
indicates T c = 38.5 K at zero magnetic field. The T c is 37.6 K for \1qH (=7.5T) || a axis 
and 36.4 K for fi H (=7.5T) || c axis. The 1/Ti was determined from an excellent fit- 
ting to 1 - M(t)/M(oa) = 0.1exp(-t/Ti) + 0.9 earo(-6t/Ti), where M(t) is the nuclear 
magnetization at time t after the saturation pulse 23]. 

Figure [H (a) shows the 75 As-NMR spectra by scanning the magnetic field at a fixed 
frequency, uo/2n = 55.1 MHz. The nuclear quadrupole frequency vq is found to be 5.1 MHz 
at 100 K which is smaller than that in the Sn-flux-grown sample (5.9 MHz) 20j. Since doping 
of K increases vq [20|, this suggests that the Sn-flux-grown crystal had a higher doping rate. 
The Knight shift K was obtained from the central transition peak and determined with 
respect to W0/7 with the nuclear gyromagnetic ratio 7 = 7.2919 MHz/T. Below T c , K is 
obtained by scanning uq at a fixed field to avoid the vortex pinning effect. Above T c , we 
confirmed that the results obtained by scanning field and scanning frequency agree well. 
The effect of the nuclear quadrupole interaction was taken into account in extracting K a . 
As shown in Fig. [1] (b), both K a and K c show a sharp decrease below T c , which indicates 




FIG. 2: (Color online) The T dependence of 1/T±. The error is within the size of the symbols. The 
dashed line shows the T 3 variation. The curves below T c are fits to a two-gap model using the 
same parameters for both directions. The inset shows the semilog plot of 1/T±T vs T c /T, which 
evidences an activation type T dependence of the relaxation.. 

spin-singlet pairing ^|. 

The main panel of Fig. 2 shows the T dependence of 1/T"i that decreases rapidly below 
T c , with the reduction over about five decades. The decrease at low T is much faster than T 3 
that is expected for a c?-wave gap. As in other materials, 1/T\ shows a "knee" shape around 
half T c , which indicates multiple gaps 4j, |5j. To see the low-T behavior more clearly, we 
plot 1/TiT as a function of inverse reduced-temperature T c /T in the inset. As can be seen 
there, \jT{T shows a very good exponential behavior below 17 K. This is strong evidence 
for a fully opened gap. 

Using the s ± -wave model and introducing the impurity scattering rate i] in the energy 

I L 

spectrum, E = to + irj [24| , but neglecting the quasiparticle damping effect for simplicity, 
we can fit the data quite well. For a sign-reversing two-gap model, as seen in Fig. 2, we 
obtain Af = 5.63 k B T c , = 1.11 k B T c , N x : N 2 = 0.85 : 0.15, where N 4 is the density 
of state (DOS) on band i, and 77 = 0.044 k B T c . For a model of three bands corresponding 
to ARPES [7], we obtain A^ = 4.7 k B T c , At = 0.96 k B T c , A3 = 4.7 k B T c , N t : N 2 : 
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FIG. 3: (Color online) T-dependence of the 75 (l/TiT). The arrow indicates T, 



N 3 = 0.44 : 0.12 : 0.44, and rj = 0.022 k B T c . The r] is much smaller than r] = 0.15 k B T c in 
LaFeAsO0.92F0.08 M an d V = 0-22 k B T c in the Sn-flux grown Ba .72K .28Fe2As2 20], meaning 
;hat the present sample is much cleaner, as supported by a small resistivity of 26/iQ at T c 
22| and the much sharper spectrum-width that is only half the value for the Sn-flux grown 
crystal. The cleanness of the present crystal is the reason for the exponential behavior of 
1/Xi at low T; impurity scatterin g br ings about finite DOS that results in seemingly power- 
law T-dependence of l/Ti 



24j-|27|. It should be emphasized that the coherence peak is 



not seen even in such clean sample, which seems hard to be explained by a s ++ -wave gap. 

Next we move to the normal state. Figure |3] shows the T dependence of 1/T]T which 
increases with decreasing T down to T c , indicating strong AF SF. The \jT{F stems from 
the magnetic susceptibility at all wave vectors. When there exists strong AF SF, one may 
assume l/T x T = {l/T l r t) AS + (1/TxT) , where (l/T^T)^ is due to the susceptibility at 
the AF wave vector Q, and (I/T1T) is due to s-band electrons and the orbital hyperfine 
interaction. Note that above T = 250K, \jT\T becomes a constant. By taking the averaged 
value of 1/TxT at T > 250K as (I/T1T) , (l/TiT) AF is then obtained. Figure 1 shows the 
ratio of the relaxation due to AF SF, (Ti) AF /(Ti) AF , which is about 2. This result indicates 
that the SF is anisotropic in the spin space, as elaborated below. 



Genera 
follows 



23| : 



ly, l/Ti is related to the transverse fluctuating internal magnetic field, 5H, as 



Y) = \fZodtcos(Lo t) (5H x (t)SH x (0) + 5H y (t)5H y (0)) 



(1) 



where (• ■ ■) denotes the statistical average. 5H is related to the fluctuating moment S of 
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Fe as <5H = A • S, where A is the hyperfine coupling tensor between the As nucleus and Fe 
spins. 



For Q = (tt, 0) and Q = (0, vr) AF SF, one has |29| 



A(vr,0) 



A 

A 



, and A(0,7r) 




A 
A 



(2) 



respectively. One therefore obtains 



AF 
b 



dtcos(u; t) (S a (t)S a (0) + S b (t)S b (0) + S c (t)S c (Q)) 



4 J-oo 
2 

^rA 2 I " dtcos(u t) (S c {t)S c {0)) 



(3) 
(4) 



Since (Sj(t)Sj(0)) (j = a,b,c) can be expressed in terms of the imaginary part of the 
susceptibility x'j through the fluctuation-dissipation theorem, dtcos(u t) (Sj(t)Sj(0)) = 



2k B T X 'j(u>o) 



the anisotropy of the relaxation can be expressed as 
Raf 



(1/T0AF 2 X »(u ,Q) + 2' 



(5) 



If x'a( u 0y Q) = Xb( u o, Q) = Xc( u o, Q)i namely, if the SF is isotropic in the spin space, then 



RTf = 1-5- 



(6) 



The observed Raf shown in Fig. H]is much larger than 1.5, which follows from Eq. (5) that 
XabiyoiQ) is larger than x'ci^OiQ) by about 50%. 

We propose that the anisotropy in the SF stems from spin-orbit coupling (SOC) that 
mixes spin and orbital freedoms so that the magnetic susceptibility bears some orbital char- 
acter, which is anisotropic. We use a two-band model [30] involving spin-orbit coupled d xz 
and d yz and calculate the anisotropy. Our theoretical study starts from a two-dimensional 
Hamiltonian: 



H = Hxb + H 



ILS- 



(7) 



Here Htb = Ek^a £ fiud^^d^cr is the tight-binding Hamiltonian with orbit index v - 
l(xz), 2(yz) and spin index a = +l(t), — 1(1)- The SOC is described by H L $ = 2X LS J2i L/ 
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FIG. 4: (Color online) T-dependence of the anisotropy of T\ due to AF spin fluctuation. The 
dashed line marks the value for isotropic AF spin fluctuation. 

The dynamical magnetic response is calculated at the random-phase-approximation level, 
with the on-site intra-orbit Hubbard interaction Hy = J2ieLattice,n Urii^rii^. The longitudi- 
nal (transverse) susceptibility Xc (x+-) a ^ Q = Q = (^O) or (0) 7r ) is 

1 

Xc(+-)(q,*'0 = z zz , xo, c (+-)(q,^n), (8) 

1 - xo, C (+-)(q, w«Ji (qj 

where \ is a 2 x 2 matrix in the orbital space with the matrix elements defined by X^ v (q, r) = 
(T T S*(-q,r)S*(q,0)> and xT_(q,r) = |<T T 5+(-q, r)£-(q, 0)). T(q) = diag(f/, 17) is the 
Hubbard interaction vertex in the spin-spin channel. The bare susceptibility Xo(<l,iv n ) can 
be easily obtained by diagonalizing the Hamiltonian shown in Eq. (7). We choose the 



3l| 



nearest-neighbor hopping integral t\ = —0.1051 eV 

The calculated magnetic anisotropy above T c with A=0.2|ti| (U = 9.7|ti|) and A=0.1|ii| 
(U = 9.87|ti|), leading to Raf ~2 at T = T c , is respectively shown in Fig. 5, which is 
in qualitative agreement with the experimental finding. Here T c =0.065|ti| is obtained by a 
self-consistent calculation of a mean-field BCS model with a gap A = 0.309|ti|. 

A particular feature we find experimentally is that the T\ ratio decreases below T c and 
it approaches the characteristic value 1.5 for the isotropic SF. Below T c , it is less trivial to 
subtract the contribution of (l/TiT) , so we simply plot the raw data as shown in Fig. [61 



We emphasize, however, that this approximation does not affect our conclusion [32j, since 
the contribution from (1/TiT) to the observed \jT{T is only 15% for H \\ a and 20% for 
H || c at T — T c . 

The asymptotic value 1.5 for (Ti) c /(Ti)a implies that the SOC effect vanishes at T — > 0. 
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FIG. 5: (Color online) Calculated magnetic anisotropy above T c due to SOC. The parameters and 
fcgT are in units of \t\\, where t\ is the nearest-neighbor hopping integral. 
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FIG. 6: (Color online) T-dependence of the T\ anisotropy below T c . The dashed straight line 
indicates the value for isotropic AF SF. 

This can happen only when the gaps are fully opened in the case of spin-singlet pairing 
with A > Als- When there are nodes in the gap function, 1/Tj at low T is governed by the 
nodal quasiparticles that are spin-orbit coupled, thereby (T 1 ) c /(T 1 ) a should resume its value 
of 2.0 at T = T c . Note also that (Ti) c /(Xi) a should become ~0.75 if the AF SF completely 
vanishes [291]. Thus, our finding of the decrease of (Ti) c /(Ti) a to 1.5 at T — > is another 
strong evidence for nodeless gap and implies that the AF SF persists in the superconducting 
state. 

In conclusion, from the NMR measurements on a clean single crystal Bao.6sKo.32Fe2As2, 
we find the long-sought exponential decrease of 1/Ti at low T, which evidences a fully- 
opened gap. In the normal state, the AF SF is anisotropic in the spin space. However, the 
anisotropy diminishes below T c and vanishes at the zero-T limit, which is a feature indicating 



S 



nodeless gap. 
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